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We study the spin-resolved transport through magnetic nanostructures in monolayer and bilayer 
graphene. We take into account both the orbital effect of the inhomogeneous perpendicular magnetic 
field as well as the in-plane spin splitting due to the Zeeman interaction and to the exchange coupling 
possibly induced by the proximity of a ferromagnetic insulator. We find that a single barrier exhibits 
a wavevector-dependent spin filtering effect at energies close to the transmission threshold. This 
effect is significantly enhanced in a resonant double barrier configuration, where the spin polarization 
of the outgoing current can be increased up to 100% by increasing the distance between the barriers. 

PACS numbers: 72.25.-b, 85.75.-d, 73.21.-b, 73.63.-b, 75.70.Ak 



I. INTRODUCTION 

The electronic properties of grapheneii^ have attracted 
in the last four years huge experimental as well as theoret- 
ical attention/^ Besides its fundamental interest as a new 
type of two-dimensional electron liquid, graphene is in 
fact regarded as a promising material for future nanoelec- 
tronic devices?^ in particular in the field of spintronicS)^ 
due to its small intrinsic spin-orbiij^ii and hyperfine 
interactions.— Indeed several recent experiment9^ i^°'^^i^^ 
have by now demonstrated spin injection and detection 
in a single layer of graphene sandwiched between fer- 
romagnetic metal electrodes and observed coherent spin 
transport over micrometer scale distances. 

Motivated by these developments, in this paper 
we focus on the problem of spin resolved trans- 
port through magnetic nanostructures in graphene. 
The studies of graphenc's transport properties 
through magnetic barriers and more complex 
structuresHa4i5^^^i9^^^^^^^ address 

the problem of controlling the confinement and the 
transport of charge carriers by means of appropriate 
configurations of an external magnetic field inhomoge- 
neous on sub-micron scales. Different types of magnetic 
nanostructures have been envisioned and their properties 
investigated, e.g., barrier a^^'^°i^^'^^ , dotsi^, wiresi^ and 
superlatticesj^ 

In all these works the spin degree of freedom has 
been completely neglected. This is justified because of 
the smallness of both Zeeman splitting and spin-orbit 
coupling^i^ii in graphene. However, the continuous and 
rapid improvements in sample preparation and experi- 
mental technology and resolution in graphene research 
call for a refinement of the theoretical analysis to incor- 
porate such finer effects. More importantly, it has re- 
cently been argued that local ferromagnetic correlations 
can be induced in graphene by several different mecha- 
nisms, e.g., proximity of a ferromagnetic insulator f22i22ii^ 
Coulomb interactions,— presence of defects,— applica- 
tion of an electric field in the transverse direction in 
nanoribbonsi^ The ferromagnetism leads to a spin split- 



ting effectively similar to a Zeeman interaction but 
of much larger magnitude. For example, a ferromag- 
netic insulator deposited on top of a graphene layer 
has been predicted to produce a spin splitting of up 
to 5 meV<^ This is comparable with the orbital en- 
ergy in a field of 1 T, which is of order of 25 meV, 
and thus may have important effects. The spin trans- 
port through ferromagnetic graphene has already re- 
ceived some attention^i^i^i2ii2^i2^ but to the best of 
our knowledge all works so far have focussed only on spin 
effects and did not consider orbital effects. Aim of this 
paper then is to fill the gap and investigate this problem 
by fully taking into account the quasiparticle's spin dy- 
namics as well as its orbital motion in an inhomogeneous 
magnetic field. 

In this paper we shall focus on the effects of Zeeman 
and exchange spin splitting and postpone the interesting 
problem of spin-orbit coupling to a future work. We show 
that the spin-resolved transmissions exhibit a strong de- 
pendence on the incidence angle, which allows in princi- 
ple for a selective transmission of spin-up and spin-down 
electrons. This effect can be qualitatively understood by 
a simple semiclassical argument. The bending of the elec- 
tron trajectories under the barrier depends on the energy 
and thus, in the presence of spin splitting, is different for 
the two spin projections. The magnitude of the polariza- 
tion that can be achieved depends on the spin splitting 
and can be very large in the presence of a large split- 
ting, as, e.g., that originating from a proximity-induced 
exchange field. Moreover, in a resonant double barrier 
configuration, the polarization can be enhanced by in- 
creasing the distance between the barriers. In this case, 
in fact, even for relatively small splitting there exists an 
energy range where the polarization reaches values close 
to one for large enough distance. 

While there exist several experimental techniques to 
produce magnetic barriers, for concreteness we shall have 
in mind the case of the fringe fields created by a ferro- 
magnetic stripe deposited on top of the graphene sam- 
ple. The magnetic field generated by the stripe is known 
analyticallyi^ The corresponding vector potential in the 
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Landau gauge can be written as A = A{x, z)y, where 



A{x,z) 



Bz{x',0)dx'— / Bx{x,z')dz' 
Jo 



, (1) 



and we assume that the system is translationally invari- 
ant in the y direction. We shall neglect possible corruga- 
tions of graphene and assume that it lies flat in the z ~ 
plane. Then the orbital dynamics is only affected by the 
perpendicular component Bz , while both Bx and Bz con- 
tribute to the Zeeman interaction. In addition we shall 
also include in the Zeeman term the effects of a possible 
exchange ficldj^i^ 

The transport properties of such structures can be cal- 
culated by means of the Landaucr-Biittikcr formalism 
in terms of the spin-resolved transmission matrix t^^i, 
which gives the probability amplitude for a quasiparti- 
cle incident on the magnetic structure from the left with 
spin projection a' to be transmitted with spin q-| 38i39,40^4i 
The matrix t^^i depends on particle's energy E and inci- 
dence angle (p (see below) and must satisfy certain general 
symmetry requirements We shall be interested in the 
spin-resolved conductances 
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and, with total conductance G = X^o-o-' Gaa' and spin 
quantization axis along the x direction, the polarization 
vector of transmitted current for a spin-unpolarizcd inci- 
dent cu 
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and their integrated value 
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The rest of this paper is organized as follows. In Sec.HIl 
we discuss the spin-dependent transmission and trans- 
port through rectangular and double resonant barriers in 
single-layer graphene. In Sec. IIIII we address the same 
problem in bilayer graphene and finally in Sec. IIVI we 
draw our conclusions. 



II. SINGLE-LAYER GRAPHENE 

Let us then start focussing on single-layer graphene. 
We shall neglect disorder and interaction effects and focus 
on a single K point, where ballistic motion of charge 
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FIG. 1: Profiles of magnetic field and exchange field consid- 
ered in this paper. In case (b) the exchange field is included 
in Bx- 



carriers in an external magnetic field B = BxX + BzZ, is 
described by the Dirac-Weyl (DW) Hamiltonian 



H = VpT 



Hf,. 



(7) 
(8) 



where p = PxX+pyij, vf ~ 10^ m/s is the Fermi velocity, 
fiB the Bohr magneton and Qs ^ 2 the effective Lande 
factor. Hspin is the sum of the Zeeman interaction due to 
the magnetic field and, possibly, the proximity-induced 
exchange splitting, with an estimated valued! h^x ~ 5 
meV (corresponding to a Zeeman interaction with a field 
of about 86 T). The vector of Pauli matrices r — T^x+Tyy 
(resp. (T = UxX+cTyy+azz) acts in sublattice space (resp. 
spin space) , and the wavefunction 5* is a four component 
object, = {'^A^,^B^,^Ai,'^Bi) (the superscript T 
denotes transposition). A = A(a;, 2;)y is the vector po- 
tential in the Landau gauge, with A(x^ z) given in Eq. 
P]). Since the Hamiltonian is translationally invariant 
in the y direction and Py is conserved, the wavefunction 
can be written as ^ = e^'^y'^ip and the DW equation re- 



in the following 



duces to a one-dimensional problem, 
we use rescaled quantities: x — s- x£b, '^x,y 
A ^ ABiB and E E^f, where B is a typical value 
of magnetic field in the problem. The magnetic length 

1/2 

is £b ^ {j§) w 26 /^yB[T] nm and the orbital energy 
scale is E^ = ^ 25\/S[T] meV, where B[T] is the 
magnetic field strength expressed in Tesla. 



A. Rectangular barrier 

First we discuss the tunneling of a spinful DW quasi- 
particlc through a rectangular magnetic barrier of width 
W with an effective Zeeman field 



B^f'f xe{w/2 



(9) 



see Fig. [ija). The value of A ranges from 0.002 for a Zee- 
man field of 1 T up to A « 0.2 if an exchange splitting 
of 5 meV is included. With spin quantization axis along 
z, the Hamiltonian ([7|) is diagonal in spin space and the 
two spin components can be treated separately. Accord- 
ingly, the wavefunctions are just two-component spinors. 



The solution to Eq. ([7]) with energy E (assumed posi- 
tive for definiteness) and spin projection a describing a 
scattering state incoming from the left can be written in 
the left and right "leads" (i.e. the non- magnetic regions 
X < -W/2 and x > W/2) as 



ij^ix < ~W/2) : 

ij^ix > w/2) 
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where ra-a-' (resp. io-o-') is the probability amplitude for a 
quasiparticle incident from the left with spin projection 
a' to be reflected (resp. transmitted) with spin a. The 
matrix Wo is given by 
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with kr 



(X) = ^E^ - [ky + A{x)] 

±k^(x)+i[ky+A{x)] 

E 

zation: 



2 and ±e±^'^(^) = 
In the leads we use the parameteri- 



ky = E sm(j)i = Esm(f>f — $, (13) 
kl = kAx < -W/2) ^ E cos (14) 
kl = k4x> W/2) ^ E cos (l}f, (15) 

where $ = A{x > W/2) is the total perpendicular mag- 
netic flux through the barrier per unit length in the y 
direction. For <I> > the emergence angle </)/ is larger 
than the incidence angle 0;. Thus a finite transmission 
is only possible if <f>i is smaller than the critical angle 
<j>c — arcsin(l — ^/E). For energy smaller than the 
threshold value Eth = \^\/2 the transmission vanishes 
for any incidence angle . ^"^'^^ This condition has a simple 
geometrical interpretation. In momentum space the dis- 
persion cone after the barrier is shifted with respect to 
the cone in the region before the barrier by — <i> along 
the ky axis. Thus if the radius of the fixed energy circle 
(the Fermi line) is smaller than $/2, the equal energy 
sections of the two cones do not overlap, implying that 
kl is imaginary for any ky, i.e., the barrier is perfectly 
reflecting. Note that neither the emergence angle 0/ nor 
the critical angle (j)c depend on spin. In the barrier region 
< 14^/2 the spin-cr wavefunction can be written asi^ 



W%{x) 



DpA-q) 



(16) 



,(17) 



where Dp{z) is the parabolic cylinder function,— q = 

I— e'^ 
V2 \A(x) + ky], Pa = ^ — I, Ec — E — aX, and a and 

b complex amplitudes. Continuity of the wavefunction 

then leads to the matching condition across the magnetic 

barrier: 



1 



cos (j)f 



1/2 
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(18) 









FIG. 2: (a) Angular plot of the transmission through a rect- 
angular magnetic barrier of width = 2 (in units of ^s) with 
A = 0.2 for spin-up (solid lines) and spin-down (dashed lines) 
quasiparticles at energy E = 1.05, 1.1, 1.7 (from left to right, 
in units of E,-n); (b) the same for a barrier of width W — 6 
at energy E = 3.05, 3.1, 3.7. Here we set B = 1 T, so that 
Em ~ 25 meV and an energy increment of 0.1 corresponds to 
2.5 meV. 



where the (spin dependent) transfer matrix T is given by 
f = [>Vo(-^F/2)]-^ WU-W/2) [WUW/2)]-' Wo{W/2). 

The transmission coefficient can directly be read from 
Eq. (UHl): 



\taa\ 



2 COS (j)f 1 



COS( 



(19) 
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The spin-resolved transmission as a function of the in- 
cidence angle is illustrated in Fig. [5] for two different bar- 
rier widths. Remarkably we find that, even though the 
critical angle (j)c is the same for both spin projections, 
there is a large difference between the spin-up and spin- 
down transmissions within an energy range A from the 
common threshold energy Eth, see Fig. [5Ja). This is due 
to the very sharp angular and energy dependences of the 
transmission onset. However, this effect becomes very 
small as soon as X/Eth ^ 1, as one can see in Fig. m^b). 
where Eth = 3. 

Before moving to more complex structures, it is in- 
teresting to briefiy consider the limit B oo, W ^ 
with fixed BW = <&, where the magnetic field profile re- 
duces to a (5 function, B{x) = ^6{x). Using the asymp- 
totic behavior of the parabolic cylinder function Dp{z) 
for p ^ oo and z ^ with z^/p finite)^ 
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FIG. 3: Illustration of the kinematics of the transmission 
through a resonant double magnetic barrier. The solid cir- 
cles represent the constant energy contours (Fermi lines) be- 
fore (resp. after) the structure (L, resp. R) or in between 
the barriers (C). In the leads L and R the circles are dou- 
bly spin degenerate. In C the vector potential horizontally 
shifts the cones along the ky axis, so that they are centered 
at kx = 0, ky = — $1. The spin splitting in C vertically 
shifts the cones upwards (resp. downwards) for spin up (resp. 
spin down). The horizontal dashed line represents the fixed y 
component of the momentum, conserved across the barriers. 
$ = $1 + $2 is the total flux through the structure (per unit 
length in the y direction). It is clear from the picture that if 
\E — \x\ < \ky + ^\ < \E + \x \ a spin-up particle propagates 
in C via an evanescent wave. 



after some lengthy algebra we obtain the compact result 
Wb{-W/2) [WB{W/2)y^ e^'xA.r.^ (20) 

where the dimensionless Zeeman coupling is now given 
by Ax = ^f^. An elementary calculation then obtains 
the transmission as 



cos^ Xz cos^ "h^jtl 



cos (jji cos 4>f 



• 2 \ 9 <t>i—4>f ' 

sm Xz cos^ ZL_iiL 



which, in contrast to the case of a barrier of finite width, 
is spin independent. We note that if one solves the scat- 
tering problem by considering the Hamiltonian ([7]) di- 
rectly with 'B'^ff ~ BW5{x)z and imposing the match- 
ing condition obtained by integrating ([7]) across the ori- 
gin with the prescription 5{x)9{x) = S{x)/2, one obtains 
Eq. (pij) with Xz replaced by Xz = 2arctan^. The 
precise functional dependence of the transmission on the 
(5-function strength depends in fact on the regularization, 
similarly to the case of an electrostatic barrier.— Since 
the use of the prescription 6{x)6{x) = S{x)/2 in the DW 
first-order differential equation has been criticized^ in 



the rest of the paper we shall use Eq. (pOI 



B. The resonant double barrier 

We now discuss the case of a resonant structure con- 
sisting of two rectangular magnetic barriers with oppo- 
site signs of the magnetic field and non- vanishing in-plane 





FIG. 4: Angular plot of the transmission for the profile in 
Fig. [Ub) with d — 6, Xx = 0.2, Xz — 0. Solid lines are for 
spin-up and dashed lines for spin-down, (a) Barriers of equal 
width Wi = W2^2 and energy E = 1.02, 1.18, 1.5 (from left 
to right); (b) the same as in (a) but for barriers of different 
width {Wi = 1 and 1^2 = 2). 



Zeeman splitting in between, as illustrated in Fig.[Tl^b): 



B 



Bz, < a; < 1^1 

B^x, Wi < X <Wi+d 

-Bz, Wi<x-d<Wi + W2 
0, otherwise 



(22) 



This profile should qualitatively model the realistic con- 
figuration of the stray field produced by a ferromagnetic 
stripe, which, in addition to the normal component, also 
contains an in-plane component B^- Inclusion of this 
component is crucial for the proper treatment of the spin 
dynamics 1^ 

First, we neglect the Zeeman term under the barri- 
ers, so that the problem is again diagonal in spin, with 
spin quantization axis along the x direction. By way of 
a simple geometric argument, illustrated in Fig. [31 we 
argue that this structure exhibits a strong wavevector- 
dependent spin filtering effect. Indeed, in the region C 
between the barriers the dispersion cones for spin-up and 
spin-down particles are equally shifted by — $i = —BWi 
along the ky axis with respect to the cones in the left (L) 
lead. Moreover the cones are also shifted (say for B^ > 0) 
upwards (resp. downwards) by the in-plane Zeeman split- 
ting, so that the radius (the Fermi momentum) increases 
(resp. decreases) by = fisBx/Em. As a result there 
exists a range of incidence angles in which the spin-down 
modes propagate via travelling waves in the central re- 
gion, whereas the spin-up modes only exist as evanescent 
waves and their transmission through the structure is 
exponentially suppressed with the distance between the 



5 



barriers. Formally, this can easily be seen from the ex- 
pression of the X component of the momentum in the re- 
gion C, kxcr = \/{E — (J^x)^ — (ky + which is real 
for a =l~ — and pure imaginary for a ~'\~ + as long as 
\E-Xx\ < \ky+'i>\ < \E+\x\. U\E-Xx\ > $/2 the trans- 
mission |t||P at any incidence angle is fully suppressed 
for large enough d and spin-up modes do not significantly 
contribute to the transport through the structure. 

The exact calculation of the transmission coefficients 
follows the same lines as for the single rectangular barrier 
discussed in the previous section. The results are illus- 
trated in Fig. 21 Fig. SJa) shows indeed that in a certain 
range of incidence angles the transmission for spin-up 
particles practically vanishes. By changing the ratio of 
barrier widths W1/W2 one can also control the position 
of the center of this interval, as illustrated in Fig. [Hb). 
The width of this range clearly depends on ■ It is then 
crucial to have a large in-plane spin splitting if one is to 
observe this effect. 

A simple closed formula for the transmission is easily 
obtained in the limit of S barriers of equal and opposite 
strengths $1 = -$2 = $: 



(cos kxad)' 



1 



(23) 




FIG. 5: (Color online) Plot of the x component of the polar- 
ization vector, Px, and, in the inset, the spin-resolved con- 
ductances Gji and Gn as functions of E, for a structure 
consisting of two 5 barriers with opposite signs of the mag- 
netic field at distance d = 5 (solid red line), 10 (long-dashed 
green line), 20 (short-dashed blue line), in units of £b- The 
Zeeman couplings are Xx = —0.0025 and Az — 0.005. We take 
B = 1 T and W = 2£b « 50 nm, and the energy is measured 
in units of Em ~ 25 meV. 



which explicitly exhibits the features discussed above and 
reproduces quite well the transmission for the case of dou- 
ble rectangular barriers. Resonances occur at k^^d = 7rn, 
with n a positive integer, where \tcra\'^ = 1. Upon increas- 
ing d, the number of resonances increases and they also 
become narrower. Interestingly, the positions of the reso- 
nances arc different for spin-up and spin-down electrons. 

Next we consider the general situation where we do not 
neglect the Zeeman splitting under the barriers. Then 
with spin quantization axis along the x direction, spin- 
flips can take place at the barriers, hence, in contrast 
to the previous case, an incident particle upon transmis- 
sion or reflection can change its spin state and the spin- 
filtering effect be spoiled. As we will see, however, the 
effect still survives close to some thresholds. 

In this case the wavefunction must be written as 



E 

(T = ± 



(24) 



where |ct = ± =t-i) are the eigenstates of ax with eigen- 
value fT = ±1 and ijj^ are two-component sublattice 
spinors. The solution to Eq. ([7]) with energy E for a 
state incoming from the left with spin projection a can 
be written in the left and right leads as 



ipix < 0) 



ipix > d) 



\/kx{x) 



yjkx{x) 



W(x) 



W(x 




(25) 
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FIG. 6: Maximum of the total polarization \P{Em)\ as a 
function of d for the same structure and parameters as in 
Fig. [5] In the inset: plot of the energy Em where the peak of 
|P(-B)| is located as function of d. 



where (5„ is the Kronecker delta and we introduce the 



(26) 4x4 matrix W given by 
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gi0T(x)gifca;f (21)2; 
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^ik^l{x)x ^-ik^i{x)x 

i<Plix) gikj;i{x)x _g—i<t>i{x}^-ik^i {x)x 



(27) 
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FIG. 7: (Color online) Plot of the spin-resolved conductances 
Gil ^^'^ ^11 fo"^ til® same structure as in Fig. [S] Here the 
Zeeman couplings are Ax — —0.2 and = 0.005. Same line- 
styles both for spin-up and spin-down. In the inset the plot 
of G|| = G11 for the same parameters. 
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FIG. 8: (Color online) Plot of the three components of the 
polarization vector, in the main figure, Py and Pz in the 
insets, for the same structure and parameters as in Fig. [7] 



with 

kxa = ^{E-<jX,)^-{ky+Ay, (28) 

_^g±#<, ^ ±k^a + ijky + A) ^ 
E - (tXx 

The transfer matrix is then given by 

where the 4x4 matrix ri(±) = e^*'^-'^^°'^ (non-diagonal 
in spin space) implements the matching conditions at the 
S barriers, and = xiO"*". For an incident particle with 
spin (T the continuity condition implies 

V / \ / 

from which the transmission amplitudes to-o-' are easily 
obtained. Then using Eqs. (HI)-® we can calculate the 
spin-resolved conductance and the spin polarization of 
the outgoing current for an unpolarized incoming cur- 
rent. The results are illustrated in Figs. OHl The 
X component of the polarization vector is plotted in 
Fig. O as a function of energy for different values of the 



distance d between the two magnetic barriers. Just for 
convenience, here and below we take B^^^ negative, so 
that the spin-up current is favored against the spin-down 
current and Px is mostly positive. For the Zeeman cou- 
plings we use the values A^; = f-LsBx/Em ~ —0.0025 and 
= \,^^^ 0.005 corresponding to -B^ = B = 1 T 
and = 2^B « 50 nm. The spin-resolved conductances 
and are plotted in the inset of Fig. [51 The con- 
ductance Gil is slightly lower than G^ while Gn = G|| 
arc negligibly small. There is then a very narrow en- 
ergy region, close io E = Eth {Eth = 1 for i? = 1 T 
and W = 2£b), in which G^ and Gn are different (see 
the inset in Fig. O, and Px exhibits a narrow peak. We 
checked that, with these values of the parameters, Py and 
Pz are of order 10""^ — 10~^, thus negligible, and Px gives 
the most important contribution to the total polarization 

\P\ = ^Jpf^pf+F^. We have calculated |P| for dif- 
ferent values of the distance between the barriers d and 
found that the polarization maximum increases with d. 
This behavior is clearly seen in Fig. [6l where the max- 
imum of \P\ is plotted as a function of d. The energy 
Em at which the polarization reaches the maximum is 
plotted in the inset also as a function of d. We observe 
that the height of the polarization peak increases with d 
and its position meanwhile shifts towards lower energy. 
From the numerical curves we can extract the following 
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behavior for \P{Em)\ and Em- 

\PiE.)\ ^ (31) 

Em =^ Eth + ^, (32) 

where Eth = 1, Ci ~ 1090, C2 ~ 1.3, a ~ 4.0 and 
/? ~ 1.3. In Fig. [S]the dots represent the exact numeri- 
cal results while the dashed line is the fitting curve. Upon 
increasing d Em approaches the transmission threshold 
relative to the first barrier. From our numerical results 
we also observe that Ci grows by decreasing |Aa;| ap- 
proximately as Ci oc lAa,!"^, while the other parameters 
in Eqs. (|3ip . ([5^ only weakly depend on the Zeeman 
couplings. For small values of \Xx\ we do not have an 
efficient spin filter since the polarization peak occurs in a 
very narrow range of energy where Gf f and are both 
very small. 

In the presence of an exchange field, instead, the ef- 
fective Zeeman interaction also includes the exchange 
contribution and it is thus much larger, A^; — —0.2. 
In this case the energy range where we get polariza- 
tion effects is widened and in the spin-resolved conduc- 
tance plot, Fig. [71 we can now clearly distinguish G|| 
from Gil- In particular, within a range of approxi- 
mately 10 meV we can have transmission of particles 
with spin-up and almost perfect reflection of particles 
with spin-down, realizing a very efficient spin filter. In- 
deed, in Fig. [5] we can see that, with a ferromagnetic 
region of width d ^ 2Q£b ^ 500 nm, particles of energy 
^/2 + Xx < E < — Xx (we set A^^ negative), i.e. be- 
tween approximately 20 and 30 meV, get perfectly spin- 
filtered upon crossing the magnetic structure. The po- 
larization, whose largest component is Px , reaches value 
one and it is sizable even for a larger range of energy. 
The other two components of the polarization vector, Py 
and Pz , due to the spin- flip processes at the barriers and 
shown in the insets of Fig. [51 remain very small, since Xz 
is small and confined to very narrow regions under the 
two thin barriers. 



III. BILAYER GRAPHENE 

In this section we consider the spin transport prob- 
lem through magnetic structures in bilayer graphene. In 
contrast to single-layer graphene, the low-energy dynam- 
ics of charge carriers in bilayer graphene is governed by 
a quadratic Hamiltonianj ^^i^^ Yet, there are important 
differences with respect to a standard 2DEG, since the 
bilayer Hamiltonian is massless and chiral, i.e., the wave- 
functions for fixed spin projection are two-component 
spinors. The effective low-energy Hamiltonian for spinful 
bilayer graphene read o^^i"^^'^^ 

TT _ ^ f (tTx- iiTyf \ , „ ^qO^ 

^ ^2^\ir^.+l7ry)^ )+^sp^n, (33 j 
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FIG. 9: (Color online) Plot of the polarization vector in bi- 
layer graphene for a structure consisting of two 5 barriers with 
opposite signs of the magnetic field, W — 21b, at distance 
d = 5 (solid red line) and 10 (long-dashed green line), in units 
oils- The dimensionless Zeeman couplings are A^, = —0.054 
and \z — 0.108. The energy is given in units of i?;, ~ 1 meV, 
for B = 1T. 



with effective mass m w 0.054 me (me is the electron 
mass in vacuum) and n^.y = Px,y + f^^.y, and H^pin 
was defined in Eq. ([5]). The possibility of inducing an 
exchange coupling by proximity of a ferromagnetic in- 
sulator has also been discussed in bilayer graphene.— 
Interestingly, in this case the electronic band structure 
can be drastically modified and a gap may open. How- 
ever, in the simplest situation, namely the bilayer sand- 
wiched between two ferromagnetic insulators with the 
same orientation of the magnetization, the only effect 
is a spin splitting and Eq. p3)) with a large in-plane 
Zeeman coupling is indeed the correct Hamiltoniani^ 
Again, using 'J = e^'^^^tp the problem is reduced to a one- 
dimensional Schrodinger equation for a four-component 
wavefunction ■(/;. As in the previous section, all quantities 
are rescaled to be dimensionless, the only difference be- 
ing that the energy scale Em = fivp/iB is now replaced 
by Eb = h'^/2me% {Eb w 1 meV for B = 1 T). Here we 
directly focus on the most interesting case of a double 
resonant barrier configuration in the limit of S barriers of 
equal and opposite strength $ = BW, with in-plane spin 
splitting between the barriers. Several different magnetic 
field profiles have also been studied in Ref. but only 
for the spinlcss case. With spin quantization axis along 
the X direction, away from the S barriers the elementary 
solutions of the Schrodinger equation for spin projection 
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FIG. 10: (Color online) Plot of the spin-resolved conductances 
Gil fi'iid Cii in bilayer graphene for the same structure and 
parameters as in Fig. [9] Same line-styles both for spin-up and 
spin-down. In the inset the plot of G|j. and G|| for the same 
parameters. 



a read 




[±k^„+i{ky+A)f 



where 



(34) 
(35) 

(36) 
(37) 



with A, = fisB^/Eb (A, « -0.054 for = -1 T) 
and Xz ~ iibBW/Ei,Ib- It is convenient to arrange the 
wavefunction ip = ^'\\ T) + V'il i) ^nd its derivative ip' in 
a eight-component vector and to write it as 



where the 8x8 matrix 'W{x) is given by 



(38) 



W{x) 



( [/|_ \ 



c/x+ c/x_ n_ 

c^I+ U[_ V[_ j 



t/j- \J -J- _ 





V 



and A is an eight-component vector of complex amph- 
tudes. The matching conditions at the positions a; = 0, d 
of the (5 barriers (i.e., continuity of the wavefunction and 
jump of its derivative) can compactly be written as 

(^T >^ , V^X >i) ^ (0- ) = (+) (Vt , ' ' ^1 ) ^ (0^) ' 



where the matrix 11, non-diagonal in spin-space, is given 

by 



17(±) = 



/ To 

TTz To TKTx 

To 

V TKTx TTz To 



(39) 



FinaUy, the transfer matrix obtains as 

f = W(0-)-i >V(0+) >V(d-)-i J7(-) W{d+). 

The scattering state for a quasiparticle of energy > 
ky and spin projection a incident on the structure from 
the left can then be written as W{x)Lc for x < Q and 
W{x)Rcr for X > d, where 

^1 = {S-\,a,rTc,,aa,0, Si^a, ri„,a'^,Q) , (40) 

^{ti^, 0,0, b„,t^„, 0,0, b'^). (41) 

are found by solving 



The transmission amplitudes tg 
the two linear systems 



TRa, 



(42) 



for cr I J,.( This can also be easily generalized to the 
case of unequal strengths of the 8 barriers.) Then, from 
Eqs. (H])-® we can calculate the spin resolved conduc- 
tance and the polarization. 

In Fig. [S] all three components of the polarization vec- 
tor are shown for two different values of d. Looking at 
the structure of the polarization, we can distinguish two 
different behaviors at two energy scales. The first occurs 
close to i? ~ Ax ~ 0.1 and is dominated by spin-flip pro- 
cesses, as one can recognize by looking at the profile of 
the three components of the polarization which all ex- 
hibit some features. Indeed, at this energy scale, \z is 
not negligible and spin-flips may play a role. At that en- 
ergy, however the conductance is almost zero, see Fig.lTUl 
so spin-flips can hardly be detected by direct transport 
measurements, at least for this value of X^. The second 
behavior occurs close to E ^ Etn = $/2 = 1 which is the 
signature of a real spin-filter effect, as one can see from 
the spin-resolved conductances plotted in Fig. (TUl At this 
energy scale Az is negligible and, in fact, Py and P^, as 
well as G\i and Gi\, which arc due to spin-flip processes, 
practically vanish, while Px reaches its maximum. 

In the presence of a larger effective in-plane Zeeman 
coupling A^;, possibly produced by the exchange spin 
splitting, the spin-filter effect is more pronounced. In 
Fig. [TT] the polarization vector is shown for two different 
value of d, with dimensionless Zeeman couplings A^; = — 1 
and Az ~ 0.1. With this larger absolute value of Aa; the 
energy range in which the spin filtering occurs is signifi- 
cantly widened. This is also seen in the behavior of the 
spin-resolved conductances in Fig. 1121 which shows in- 
deed that particles with spin-down are almost perfectly 
reflected by the magnetic structure for energies smaller 
than approximately 2 — 3 meV, while spin-up particles 
are almost always transmitted with spin up, since the 
amplitude for spin-flips is very small, as shown in the 
inset of Fig. [H 
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FIG. 11: (Color online) Plot of the polarization in bilayer 
graphene, for a structure consisting of two S barriers with 
opposite signs of the magnetic field, W = 2iB, at distance 
d = 5 (solid red line) and 7 (long-dashed green line) , in units 
of ^s. The dimensionless Zeeman couplings are = — 1 and 
— 0.108. The energy is given in units of Et ^ 1 meV for 

B = n. 

IV. CONCLUSIONS 

In conclusion, we have analyzed the spin transport 
problem through magnetic nanostructures in graphene. 
We have shown that an inhomogeneous field profile to- 
gether with a strong in-plane spin splitting can produce 
a remarkable wavevector-dependent spin filtering efi^ect. 
This effect is enhanced in a resonant barrier configura- 
tion, where the polarization can reach values up to one. 
This result can be understood by means of a simple kine- 
matical analysis of the problem. 

While we confined ourselves to zero temperature, we 
expect that the effect should be observable at finite tem- 
perature as well, as long as the temperature is smaller 
than the in-plane spin splitting. If the splitting origi- 
nates from the exchange coupling and the estimate of 
Ref. is experimentally confirmed, then there is a com- 
fortable temperature window (say below 10 K) where 
the effect we discussed could in principle be observed. 
Other mechanisms inducing local ferromagnetic correla- 



tions in graphene could also be exploited to increase the 
spin splitting, thereby improving the spin filtering effect. 
Moreover, since the orbital dynamics and the spin dy- 
namics in this problem are to a large extent decoupled 
(this would not be the case with spin-orbit coupling), we 
expect that the addition of a small amount of impurity 



0.4 
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FIG. 12: (Color online) Plot of the spin-resolved conductances 
Gil ^11 ill bilayer graphene, for the same structure and 
parameters as in Fig. 1111 Same line-styles both for spin-up 
and spin-down. In the inset the plot of G11 and G11 for the 
same parameters. 



scattering would not spoil the spin filtering effect, at least 
as long as the scatterers are enough long-range that they 
do not induce scattering between the two K points. 

Along the same lines of this work, one could also in- 
vestigate the effects of spin-orbit coupling (SOC). While 
the SOC has so far been estimated to be very small in 
graphenef^ii recent experimental results^ indicate that 
in quasifreestanding graphene produced on Ni(lll) with 
intercalation of Au the Rashba effect leads to a large 
spin splitting of order of 13 mcV. We plan to address 
this problem in a forthcoming work. 

Finally, we hope that our paper will stimulate further 
experimental research on the physics and the transport 
properties of magnetic nanostructures in graphene. 
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